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Using ultrafast optical pump-probe spectroscopy, we study the relaxation dynamics of hot optical 
phonons in few-layer and multi-layer graphene films grown by epitaxy on silicon carbide substrates 
and by chemical vapor deposition on nickel substrates. In the first few hundred femtoseconds after 
photoexcitation, the hot carriers lose most of their energy to the generation of hot optical phonons 
which then present the main bottleneck to subsequent carrier cooling. Optical phonon cooling on 
short time scales is found to be independent of the graphene growth technique, the number of layers, 
and the type of the substrate. We find average phonon lifetimes in the 2.5-2.55 ps range. We model 
the relaxation dynamics of the coupled carrier-phonon system with rate equations and find a good 
agreement between the experimental data and the theory. The extracted optical phonon lifetimes 
agree very well with the theory based on anharmonic phonon interactions. 
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The performance of most demonstrated and proposed 
graphene-based electronic and optical devices depends 
critically on carrier and optical phonon scattering as 
well as relaxation dynamicsi^^. Optical phonons in 
graphene can decay in two main ways: i) optical phonons 
can get absorbed by electrons or holes, or ii) optical 
phonons can decay into other phonons via anharmonic 
interactions and lattice defects. In this paper, phonon 
lifetime refers to the lifetime due to the latter processes, 
the most important of which is the decay of an optical 
phonon into two acoustic phonons^. Optical phonon life- 
times in graphene can be estimated from frequency do- 
main measurements, such as Raman spectroscopy. How- 
ever, the measured Raman linewidths also contain con- 
tributions from pure dephasing processes and inhomoge- 
neous broadening 7 . Direct measurement of the optical 
phonon lifetime in carbon nanotubes via time-resolved 
Raman anti-Stokes spectroscopy was performed by Song 
et. all 7 - and a value of ~1.1 ps was reported. Theoret- 
ical calculations of the hot optical phonon lifetimes in 
graphene due to anharmonic decay into acoustic phonons 
were carried out by Bonini et. al£, and values between 
2-3 ps were reported for both zone-center (E2 g ) and zone- 
edge (A[) modes for phonon temperatures in the 500-900 
K range. A much longer lifetime value of ~7 ps was 
measured by Kampfrath et. al. 8 via time-resolved mea- 
surements of the electronic temperature using optical- 
pump terahertz-probe spectroscopy. Degenerate opti- 
cal pump-probe studies of carrier dynamics in epitaxial 
graphene have been reported previously by the authors 
and other a 12 ' 14 . In this paper, we use optical pump- 
probe spectroscopy to study optical phonon lifetimes in 
graphene, and we also present a theoretical framework 
for interpreting the results. 

Graphene samples used in this work were grown on 
the silicon and carbon faces of semi-insulating silicon 
carbide (6H-SiC) wafers by thermal decomposition (epi- 
taxial growth^) and also by chemical vapor deposition 



(CVD) on nicke l 10 ' 11 . Samples A and B were grown 
on the carbon face of SiC and Sample C was grown on 
the silicon face of SiC. Sample D was grown by CVD 
on nickel and then transferred onto a quartz substrate. 
Sample D was found to have patches in which the num- 
ber of layers varied between 1 and 3. All samples were 
characterized via Raman spectroscopy (pump wavelength 
488 nm) and optical/IR transmission spectroscopy. The 
characteristics are recorded in Table I. For spectroscopy, 
the pump and the probe beams were obtained from a 
TkSapphire mode-locked laser with a 81 MHz pulse rep- 
etition rate and a 780 nm center wavelength. The pulse 
width (FWHM) at the sample was measured to be ^100 
fs. The polarization of the probe was rotated 90 with 
respect to the pump such that scattered pump light could 
be filtered out of the probe beam before photodetection 
using a polarizer and a spatial filter. The differential 
transmission (DT) of the probe due to the pump was 
obtained by chopping the pump and probe at different 
frequencies and recording the photodetector current at 
the sum frequency with a lock-in amplifier. 

The measured time-resolved DT signals (normalized to 
their peak values) for all the four samples are shown in 
Fig.l (solid circles). In each case, the pump energy was 
~14 nJ (ffuence ~15 /iJ/cm 2 ). The photoexcited carrier 
density was estimated to be ~5xlO n cm -2 per layer 
for Samples A through C (SiC substrate), and ~8xlO n 
cm~ 2 per layer for Sample D (quartz substrate). The 
autocorrelation of the ^100 fs pulses used in the experi- 
ments are also shown in Fig.l (dashed lines) for reference. 
The measured transients exhibit two distinct time scales. 
The DT signal decreases rapidly to almost 10% of its 
peak value in the first ^400-500 fs, and then it decreases 
much more slowly afterwards. These results can be qual- 
itatively explained in the following way. We assume that 
the probe transmission is affected only by the change 
in the interband absorption caused by the pump pulse. 
Since the probe transmission is sensitive to the carrier 
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FIG. 1: Measured differential probe transmissions (normal- 
ized to peak value) (circles). The autocorrelation of the ~100 
fs pump and probe pulses (dashed lines), and the best- fit the- 
oretical curves (solid lines) are plotted for different graphene 
samples. The extracted values of the average optical phonon 
lifetime r v u are indicated for each graphene sample. 



occupation in the conduction and valence bands at en- 
ergies ~%u>/2 above the Dirac points (half of the pump 
photon energy) , it is also sensitive to the carrier temper- 
ature. Immediately after photoexcitation, the photogen- 
erated electrons and holes thermalize with each other and 
the existing carriers, thereby acquiring Fermi-Dirac dis- 
tributions with high temperatures. Recent studies have 
shown that the thermalization times for the photoexcited 
carriers in graphite are extremely short and in the 20-40 
fs ranged. Due to the fast electron-hole scattering times, 
it is also reasonable to assume that the electron and hole 
temperatures are approximately the same. The thermal- 
ized electron and hole distributions then cool down via 
interaction with the optical phonons. The hot photoex- 
cited carriers lose most of their energy to the optical 
phonons in the first 500 fs after photoexcitation. This 
results in the generation of hot optical phonons, which 
then present a bottleneck to the subsequent cooling of 
the carriers. The relaxation dynamics of the carriers and 
the optical phonons are strongly coupled, and experimen- 
tal data can be interpreted correctly only if this coupling 
is taken into account. Below, we present a model that 
considers this coupling, and we use this model to extract 
the optical phonon lifetime from the data. 

Intraband electron-phonon scattering can be intraval- 
ley (due the T-point E2 g optical phonons) or intcrvalley 
(due to the K(K')-point A' x phonons). Using the matrix 
elements for electron-phonon scattering given by Rana 
et. alii£, the net optical phonon emission rates, Rr c and 
Rkc (units: cm _2 s _1 ), due to intraband intravalley and 
intervalley electron-phonon scattering, respectively, can 



TABLE I: Experimentally determined quantities: Ratio of 
Raman G-peak and D-peak intensities Ig/Id, the number of 
graphene layers N, the extracted values of the optical phonon 
lifetime r p h and the average carrier density per layer n . 
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Tph 




n c 
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20.5 


55 


2.55±0.08 


ps 


lxlO 11 


cm" 
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57 


16 


2.5±0.08 


ps 


lxlO 11 


cm" 
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9.6 
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2.55±0.1 


ps 


6xlO n 


cm" 
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19 


~2 


2.5±0.12 


ps 


7x10" 


cm" 
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FIG. 2: (a) The electron intraband scattering time due to op- 
tical phonon scattering is plotted as a function of the electron 
energy^, (b) The phonon Brillouin zone is shown together 
with the annular regions (not to scale) at the F and K(K') 
points, in which the optical phonons participate in electron- 
phonon scattering. 



be expressed as, 
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The integrand in the above Equation is the distribution 
in energy of the electron-phonon scattering events. The 
expression for i?Ke is the same as above except for the re- 
placement r — > K. Here, dt/db equals ~45 eV/n m 16 ' 17 , 
p is the density of graphene and equals ~ 7.6 x 10~ 7 
kg/m 2 , (nr) and ojk (j*k) are the frequencies (av- 
erage occupation numbers) of the Ei g and A\ phonons, 
respectively, v = 10 6 m/s, and f c (E) is the Fermi-Dirac 
distribution function for conduction band electrons. Sim- 
ilar expressions as above can be written for the optical 
phonon emission rates, i?rh and i?Kh, due to hole-phonon 
intraband scattering. The carrier temperature T c and the 
optical phonon occupation numbers, nr and nx, obey the 
coupled rate equations, 



dT c 
dt 

dnr 
dt 

dnyi 
dt 



(Rr c + Rr%)huJr + (Rk c + -BkiQ^k 
C e + Ch 
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FIG. 3: (a) Calculated electron-phonon scattering rate versus 
electron energy is plotted for values of T c equal to 1250K, 
1500K, 1750K, and 2000K. Total electron density is assumed 
to be 7x 10 11 cm -2 . Phonon temperature is fixed at 900K. (b) 
Measured differential probe transmissions (normalized to the 
peak value) (circles) and the best-fit theoretical curves (lines) 
are plotted for different pump pulse energies (2.5, 6, 9.8 nj). 
(c) and (d) The calculated electron and phonon temperatures 
and the phonon occupation numbers are plotted as a function 
of time after photoexcitation. The photoexcited electron-hole 
densities are assumed to be 2 x 10 11 cm -2 and 7 x I0 11 cm -2 . 
The equilibrium electron density is assumed to be I0 11 cm -2 . 



Here, C e and Ch are the electron and hole heat capac- 
ities and T p h is the average lifetime of the E2 g and A' x 
optical phonons. n r and ri K are the phonon numbers 
in equilibrium (at T=300K). Mr and Mk are the num- 
ber of phonon modes (per unit area) at the T and K(K') 
points, respectively, that participate in carrier-phonon 
interaction. The values of Mr and Mk can be estimated 
as follows. At the T-point, only phonons with wavevector 
magnitudes q in a certain range (dictated by momentum 
and energy conservation) participate in carrier-phonon 
interaction. If an electron with wavevector k emits an 
optical phonon, then the smallest value of q is u>r/v and 
most of the emitted phonons have wavevectors smaller 
than k. For carrier densities in the 10 11 — 10 12 cm' 2 
range, and carrier temperatures in the 300-3000K range, 
most of the electrons that emit optical phonons have en- 
ergies less than E max « 0.8 eV (see Fig. 3(a)). Therefore, 
the largest value of q can be approximated by E max /Kv. 
These estimates have been verified by detailed calcula- 
tions using wavevector dependent distributions for the 
phononsiS,. The value of Mr then corresponds to the 
number of phonon modes that are in the annular region 
(see Fig.2(b))defined by the minimum and maximum val- 
ues of q and equals 2 x [{E max /hv) 2 — {ojt / v) 2 ] / Air . The 
factor of 2 in the front stands for the two degenerate LO 
and TO phonon bands at the T-point. Similarly, Mk 
equals 2 x [(E max /hv) 2 — (luk/v) 2 ]/4:Tt. Here, the factor 



of 2 in the front stands for the two degenerate phonon 
valleys at the K and K' points. Since the optical phonon 
emission time for an electron at an energy Uu>/2 ps 0.8 eV 
above the Dirac point is less than ~50 fs (see Fig. 2(a)), 
the initial carrier temperature is computed by assuming 
that during the first 100 fs after photoexcitation all the 
photogenerated carriers emit ^2 optical phonons. The 
set of coupled Equations @-© were solved to obtain 
the carrier temperature and the phonon numbers as a 
function of time. It is assumed that no significant carrier 
recombination occurs in the first few picoseconds^. The 
carrier distributions were then used to compute the opti- 
cal transmission using the expressions given by Dawlaty 
et. ali^ 2 - and the results were convolved with a probe pulse 
assumed to be Gaussian with a FWHM of 100 fs. Theo- 
retical results obtained are plotted in Fig.l (solid lines). 
The agreement with experimental data is seen to be very 
good, and both the fast and slow time scales in the DT 
signal are correctly reproduced. It needs to be pointed 
out here that the only two fitting parameters in the model 
are the phonon lifetime, r p h, and the equilibrium carrier 
density per layer, n Q . Their extracted values are given in 
Table I. The error bars in Table I for the extracted val- 
ues of Tph correspond to the spread in the data obtained 
when the measurements were repeated many times. The 
extracted values of t p ], are not sensitive to the extracted 
values of n Q nor to the photoexcited electron-hole densi- 
ties. Fig. 3(b) shows the agreement between the theory 
and the data for pump pulse energies of 2.5, 6, and 9.8 nJ 
(fluences 6.2, 15, and 24.5 /iJ/cm 2 ). In the model, the 
photoexcited carrier density was scaled with the pump 
pulse energy. In each case, the average best fit value for 
Tph was found to be 2.5 ps. 

Fig. 3(c) and Fig. 3(d) show the calculated carrier and 
phonon temperatures and phonon occupation numbers 
plotted as a function of time after photoexcitation for 
photoexcited electron-hole densities of 2 x 10 11 cm~ 2 and 
7 x 10 11 cm~ 2 . Tph was assumed to be 2.5 ps. We as- 
sumed that the phonons were initially at room temper- 
ature (300K) and the equilibrium electron density was 
10 11 cm~ 2 . In the first 500 fs, the electron temperature 
rapidly decreases and the phonon temperature rapidly 
increases. For photoexcited carrier densities in the 10 11 - 
2xl0 12 range, the maximum phonon temperatures were 
found to be in the 550-1550 K range, in agreement with 
the earlier findings of Song et. As the tempera- 

ture of the carriers and the phonons approach each other, 
the net energy exchange between the electrons and the 
phonons also decreases. The hot optical phonons then 
become the main bottleneck for further carrier cooling. 
It should be noted that the subsequent cooling of the 
coupled carrier-phonon system is not determined by the 
phonon lifetime alone. The coupling of the hot carriers to 
the phonons retards the cooling of the phonons. Fig. 3(d) 
shows that the phonon occupation numbers exhibit de- 
cay times that are slightly longer than the value of T p h 
and in the 2.9-3.5 ps range. Because the temperatures of 
the r and K(K') point optical phonons remain close to 
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each other, larger-energy T-point phonons exhibit slightly 
smaller occupation numbers and faster decays compared 
to the K(K')-point phonons. 

The results from Samples A through D show that the 
optical phonon cooling rates on short time scales are inde- 
pendent of the growth technique, the number of graphene 
layers, and the type of the substrate. These results imply 
that on the fast times scales relevant to our experiments 
energy transport among the graphene layers, or between 
the graphene layers and the substrate, is not the bottle- 
neck for hot optical phonon relaxation. The values of 
the optical phonon lifetime measured in this work agree 
very well with the theory based on anharmonic phonon 



interactions developed by Bonini et. ali£. Contrary to 
expectations^, the substrate (SiC and quartz) optical 
phonons do not seem to play noticeable or distinctive 
role in hot carrier relaxation. This could be because the 
samples used in this study were not single-layer. 
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